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Abstract: We prove that various SO (n) -invariant n-vector mod- 
els with interactions which have a deep and narrow enough mini- 
mum have a first-order transition in the temperature. The result 
holds in dimension two or more, and is independent on the nature 
of the low-temperature phase. 

Recently Blote, Guo and Hilhorst , extending earlier work by Domany, 
Schick and Swendsen on 2- dimensional classical XY-models, performed 
a numerical study of 2-dimensional n-vector models with non-linear inter- 
actions. For sufficiently strong values of the non-linearity, they found the 
presence of a ffist-order transition in temperature. In ^ a heuristic explana- 
tion of this ffist-order behavior, based on a similarity with the high-g Potts 
model, was suggested, explaining the numerical results. A further conffi- 
mation of this transition was found by Caracciolo and Pellisetto M, who 
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considered the n — > oo (spherical hmit) of the model, and found the same 
first-order transition. 

On the other hand, various studies, mostly based on Renormalization- 
Group analyses or Kosterlitz-Thouless type arguments based on the picture 
of binding-unbinding of vortices, have contested this first-order behaviour 
(e.g. 

Here we settle the issue by presenting a rigorous proof of the existence of 
this first-order transition. 

It may seem somewhat surprising that two-dimensional n-vector models, 
whose magnetization by the Mermin- Wagner theorem |^ is always zero, can 
have such a phase transition. The reason is that the transition we are talking 
about here is manifested by the long-range order in higher-order correlation 
functions. Such transitions were discovered by one of us some time ago, see 
|rT|. But the results of |TT| were related to the fact that there the symmetry 
group was the (disconnected) group O (2) , and at the transition point the 
discrete symmetry Z2 is broken, while the connected part - SO (2) - of the 
symmetry persists. The nature of the transition we study here, however, is 
not connected to any type of symmetry breaking, and as such is much closer 
to the high-g g-state Potts model, or the model studied in 0], where a first 
order transition between a low-energy and a high-entropy phase occurs. Thus 
we confirm the original intuition of 

For XY-spins in two dimensions there can be a low-temperature phase 
with slow polynomial decay of correlations, while the majority belief in the 
field, despite the work of Patrascioiu and Seller 0, is that for n > 2, the 
n-vector models at low finite temperatures have exponentially decaying cor- 
relations, just as at high temperatures. Our result unfortunately does not 
say anything about this question. 

Our proof is directly inspired by the existing proofs for low-energy-high- 
entropy phase transitions, and is indeed an adaptation of those. We employ 
the method of Reflection Positivity |T^ (RP). For simplicity we write the 
proof for 2-dimensional XY-spins, the extension to the general case is imme- 
diate. 

We remark that also the generalization to higher dimensions is immedi- 
ate. Thus the low-temperature phase can be either magnetized, Kosterlitz- 
Thouless-like (not magnetized with slow correlation decay), or possibly non- 
magnetized with exponentially decaying correlations. As such our result 
contradicts strong "universality" claims, stating that universality classes of 
interactions exist, all elements of which have the same kind of phase tran- 
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sition between a high-temperature phase and a low-temperature phase, and 
which are only determined by the dimension of the system, the symmetry of 
the interaction and whether the interaction is short-range or long-range. 
We will consider the Hamiltonian given by 



To formulate our result we have to introduce for every n.n. bond b = 
the following bond observables: 

n {<P^,<PJ) -jo if 10,-0,1 >£/2, 



which project on the ordered bond configurations, and 
(0j, 0j) . Our main result is contained in the following 

Theorem 1 Suppose the parameter p is large enough. Then there exists a 
transition temperature [3c = Pc {J,p) , such that there are two different Gibbs 
states, (■)^ and (■)^ , at P = Pc, corresponding to the Hamiltonian For 
some specific choice of e = e (p) in we have for the "ordered" state (■)^ 
that 

{P,^f>^{p), 
while in the "disordered" phase (■)^ 

for each bond b, with k (p) — > 1 as p — > cx). 

Before analysing the model (Q), we present an even simpler toy model, 
which already displays the mechanism, and which is even closer to the Potts 
model. The single spin space is the circle, S^, the free measure is the Lebesgue 
measure, normalized such that has measure one, so = [— |, . ( One 
can take here any sphere §" instead.) The toy Hamiltonian is 

H = -J J2 U (3) 
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The nearest neighbour interaction U (0i, ^2) = (|0i — 02|) is given by 



-1 if 101 < 



£ 

2' 

otherwise. 



Here e plays a similar role to ^ in the g-state Potts model. The Hamiltonian 
is RP under reflections in coordinate planes. For the case one has also 
RP under reflections in lines at 45 degrees, passing through the lattice sites. 
That case is the easiest. 
One has to show that: 

• {Pi^) p is small for large [3] (the ordered, typical low-temperature phase 
bonds) 

• {P^)p is small for small f3; (the disordered, typical high-temperature 
phase bonds) 

• {PyPy,)p is small for all /3, 

provided e is small enough. Here (■)^ is the state with periodic b.c. in 
the box A of size L, and 6', h" are two orthogonal bonds sharing the same 
site. The estimates have to be uniform in L, for L large. The first two are 
straightforward application of RP and the chess-board estimate. So let us 
get the last one. By the chess-board estimate, 

(n^n^), < {p*)T ' (4) 

where the observable is the indicator 



^* = n A" n 



Here i?oi,-E'23 is the partition of all the bonds in A into two halves; Eqi 
consists of all bonds (x, x + ei) and {x,x + 62) , for which + = or 
Imod 4, while i?23 is the other half; ei and 62 are the two coordinate vectors. 

To proceed with the estimate (^) we need the estimate on the partition 
function. We have 

^A(Ae)>(|)'"e^^l^l + (l-4^)l^l/^ (5) 
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(The first summand is obtained by integrating over all configurations 0, such 
that 10^1 < I for all x G A. For the second one we take all configurations (j) 
which are arbitrary on the even sublattice and which satisfy \4>x — 4>y\ > sj2 
for every pair of n.n.; for every y on the odd sublattice that leaves the spins 
to be free in a set of measure > 1 — 4£.) Solving 



for /3o we find 



4£) 



|A|/2 



-2 



(6) 



so for P > Po the first term in dominates, while for P < Po the second 
term dominates. Similarly, the partition function Zf {P, e) , taken over all 
configurations with (0) = 1 satisfies 



If P > Po, we write, using @: 



(7) 



2^ 



£4 6^ 



< 2- 



If /? < /3o, we similarly have 



\PuPy')0 — 



;i-4£) 



1/2 



< 



5(1-4£)(|) 



-2 



(l-4£) 



1/2 



So we are done. 

For the non-linear models, we employ the fact that for small difference 
angles cos(0i — 0j) is approximately 1 — O((0j — 0j)^) and furthermore that 



lim^ 



e ^. This suggests to choose e(p) = -j=. Because the 



ip^Oo(l pY — ^ ■ -^^^^^ ^^^^^^^^ ^^^^^^■^ ^y^yj — 

separation between ordered and disordered bonds is somewhat arbitrary, to 
obtain an inequality similar to (5), we make a slightly different choice. We 



consider a bond (i, j) disordered if |0j — ( 



> 



c_ 

VP 



for some large C. So first we 
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choose a sufficiently large constant C. For the estimate of the ordered parti- 
tion function we only integrate over the much smaller intervals of "strongly 



ordered" configurations: < to obtain a lower bound: 

/ 1 \|A| / 4^\|A|/2 

This makes use of the fact that the "strongly ordered" bonds all have 
energy almost equal to —J, whereas the disordered partition function is 
bounded by that of the toy model, but with e replaced by 

For the estimate which shows that ordered and disordered bonds tend 
not to neighbor each other, we obtain: 



The rest of the argument is essentially unchanged. We first chose C big 
enough (such that ^ is small wrt to 1), and we can still choose p big enough 
such that p is large wrt to C^, which finishes the argument. 
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